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Abstract. We investigate the role of rotational effects in inducing asymmetry 
present above ~ 5 GeV in the double-peak lightcurves of the bright EGRET pul- 



sars: Vela, Crab, and Geminga. According to Thompson (2001), the trailing peak 
dominates over the leading peak above ~ 5 GeV consistently for all three pul- 
sars, even though this is not the case over the entire energy range of EGRET, i.e. 
above ~ 100 MeV. We present the results of Monte Carlo simulations of electro- 
magnetic cascades in a pulsar magnetosphere within a single-polar-cap scenario 
with rotationally-induced propagation effects of the order of (3 (where (3 is the 
dimensionless local corotation velocity). We find that even in the case of nearly 
aligned rotators with spin periods of P ~ 0.1 s rotation may lead to asymmet- 
ric (with respect to the magnetic axis) magnetic photon absorption which in turn 
leads to asymmetric gamma-ray pulse profiles. The resulting features - softer spec- 
trum of the leading peak and the dominance of the trailing peak above ~ 5 GeV 
- agree qualitatively with the EGRET data of the bright gamma-ray pulsars. 
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1. Introduction 

Good quality gamma-ray data for three pulsars - Vela, Crab, and Geminga - acquired 
with EGRET aboard the CGRO tempts to analyse the properties of pulsar high-energy 
radiation as a function of photon energy and phase of rotation. Gamma-ray spectra of 
pulsed radiation from these sources (as well as from three other EGRET pulsars: B1706- 
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44, B1951+32, and B1055-52) extend up to ^ 10 GeV. All three pulsars feature gamma 
lightcurves characterised by two strong peaks separated by 0.4 to 0.5 in rotational phase. 
These double-peak pulses are asymmetrical and their profiles change with energy. Above 
~ 100 MeV the leading peak (LP) is stronger than the trailing peak (TP) in the case of 
the Vela and the Crab pulsars, and the opposite is true for Geminga. However, for all 
three pulsars their leading peaks exhibit lower energy cutoffs - around ^ 5 GeV - than the 
trailing peaks (TP). In other words, the trailing peaks dominate over the leading peaks 



above ~ 5 GeV (Thompson 2001 ) . In the case of the Vela pulsar and Geminga, this effect 



is accompanied by the softening of the spectrum of the leading peak (Fierro et al. 199g , 



Kanbach 1999 ). Potential importance of the double-peak pulse asymmetry in the case of 
Vela was already acknowledged - at the time when the COS-B data became available - 



by Morini 1983 who attempted to explain the asymmetry with a hybrid model, with two 
different mechanisms responsible for the formation of the leading and the trailing peak. 

High-energy cutoffs in pulsar spectra are interpreted within polar cap models as due 
to one-photon absorption of gamma-rays in strong magnetic field with subsequent e^- 
pair creation. A piece of observational support for such an interpretation comes from 
a strong correlation between the inferred 'spin-down' magnetic field strength and the 



position of the high-energy cutoff (Baring & Harding 2000, Baring 2001). This in turn 
opens a possibility that the observed asymmetry between LP and TP, i.e. the dominance 
of LP over TP above ~ 5 GeV, is a direct consequence of propagation effects (which 
eventually lead to stronger magnetic photon absorption for photons forming LP than 
TP) rather than due to some inherent property of the gamma-ray emission region itself. 

The aim of this paper is to investigate the role of pulsar rotation in a built-up of 
such asymmetry in the double-peak pulse profiles. We consider purely rotational effects: 
due to presence of rotation-induced electric field E, aberration of photon direction and 
slippage of magnetosphere under the photon's path. In section 2 we compare them with 
some other effects which may be responsible for the asymmetry (like various distortions 
of the magnetic field structure). In section 3 we show that the rotation effects result in 
an asymmetric pair production rate for the leading and the trailing part of the magneto- 
sphere even in the case when the magnetic field structure and the population of radiating 
particles are symmetric around the magnetic pole. In section 4 asymmetric pulse pro- 
files are calculated as a function of photon energy and then the model predictions of 
the ratio of fluxes in the leading and trailing peaks are compared with the inferred ratio 
for Vela at different energy bins. In section 5 we address the significance of rotation- 
driven asymmetry across the pulsar parameter space. Our main results are discussed in 
Section 6. 
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2. Symmetric features of the model 

It has been argued in many studies of radio properties of pulsars (eg. Blaskiewicz et 
al. 1991, Gangadhara & Gupta |2001| ) that a rigidly rotating static-like dipole can be 
used as a good approximation of dipolar magnetic field as long as only most important 
rotation effects, of the order oi j3 — v/c (where v is the local corotation velocity and c is 
the speed of light), are to be considered. According to order-of-magnitude estimates, small 
distortions of the dipolar magnetosphere induced by rotationally-driven currents can be 
neglected: longitudinal currents suspected to flow within the open field line region cannot 
modify B by a factor exceeding whereas toroidal currents due to plasma corotation 
change B barely by a value of the order of . A more comprehensive discussion of the 
influence of currents on the magnetospheric structure can be found in Beskin ( |1999|) and 
references therein. 

Below we follow the approximation of a rigidly rotating static-like dipole: at any 
instant the magnetic field has the shape of a static dipole in the frame which corotates 
with a neutron star. Moreover, the magnetosphere is assumed to be filled out everywhere 
with the Goldreich- Julian charge density, so that a rotation-induced electric field E = 
— /3 X i? is present in the OF, whereas E' = in the frame corotating with the star. We 
neglect deviations from this corotational electric field which are present in the charge- 
deficient polar gap region - they do not exceed a factor of /J'^/^ . 

In our Monte Carlo simulations (section 4) development of gamma-ray radiation is 



based on the model of Daugherty & Harding |1982[ with primary electrons being injected 
along magnetic field lines at an altitude of a few neutron star radii, at the magnetic 
colatitude corresponding to the last open magnetic field lines, and uniformly in the 
magnetic azimuth. The electrons are assumed to accelerate instantly to the energy Eq ^ 
10'' MeV and subsequently to cool down emitting curvature photons. Some of the photons 
induce in turn electromagnetic cascades which propagate outwards in a form of a hollow 
cone beam (see Dyks & Rudak 200(\ for detailed description of directionality aspects of 
the casades as well as viewing geometry). 



3. Rotation-driven asymmetry 

In this section we present in detail the mechanisms which lead to rotation-induced asym- 
metry in the (otherwise axially symmetric) hollow-cone gamma-ray beam. For the sake 
of better demonstration of the effects we consider an exaggerated case: propagation of 
a photon in the equatorial plane of a fast orthogonal rotator. Fig. |l| shows the case for 
the rotation period P = 1.5 ms: photons of the same energy in the corotating frame 
(CF) are emitted from two opposite points on the outer rim of the polar cap tangen- 
tially to the magnetic field lines in the CF. To follow their straight-line propagation in 
the inertial observer frame (OF) three effects have to be taken into account: 1) In the 
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Fig. 1. Top view of orthogonally rotating pulsar with the spin period P — 1.5 ms; the 
sense of rotation is counter-clockwise (indicated) ; the cartesian coordinates are in units of 
the light cylinder Ric = cP/2'k. The magnetic field lines are approximated with the static- 
like dipole of the magnetic moment /x (indicated). Two photon trajectories (starting from 
two opposite points on the polar cap) in the inertial observer frame (OF) are marked 
with two long dotted lines. These are the prototypes of the leading peak (on the left) and 
the trailing peak (on the right). In the corotating frame (CF) these trajectories follow 
the long dashed lines. Short segments of the dotted-line trajectories (OF) have been 
transformed via rotation to several actual photon locations in the magnetosphere, and 
thus they indicate the local propagation directions fj of the photons in OF. Solid arrows 
at these locations point along the directions i)''^ of local free propagation in OF. Note 
much larger angles between 77^^ and r) for the photons of the leading peak than for the 
photons of the trailing peak. 

OF the photons propagate at an aberrated direction (dotted lines in Fig. |l|) and differ 
in energy. 2) The rotation-induced electric field E = —(3 x J5, which is present in the 
OF, modifies the rate of the magnetic photon absorption in a different way for photons 
forming the LP than for photons forming the TP. 3) The electromagnetic field in the OF 
is time-dependent: because of the rotation the photons propagate through different parts 
of the magnetosphere. 

We find that the second effect - due to rotationally induced E - plays a dominant role 
in generating the asymmetry in the magnetic absorption rate R between photons of the 
LP and the TP. An importance of a weak electric field E J- B for the rate R was for the 
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first time recognized by Daugherty & Lerche (1975) who presented also its quantitative 
treatment. A consequence of its presence is that the rate R does not vanish along the 
direction of local B any more, and becomes non-axisymmetric with respect to it as well. 
Instead, the rate R vanishes along two new directions which lie in the plane perpendicular 
to E and deviate from B by angle ~ E/ B in such a way that in a local coordinate frame 
with z \\ B, y \\ E, and x \\ Ex B the "free propagation" direction 77^^ — [Vx'' iVy'' iVz''] 
has two solutions: [E/B, 0, ±(1 — E^ / B^y^'^]. For definiteness, hereafter we will consider 
photons which propagate outwards within the regions above the northern magnetic pole 
(i.e. with propagation vectors 77 satisfying i) ■ B > Q) which corresponds to the case of 
77^^ = +(1 — E'^ / B'^y/'^ . Fig. H shows the mean free path Amfp = c/R for the magnetic 
photon absorption for different directions in the plane perpendicular to E. The rate R was 
calculated by performing Lorentz transformation to a frame in which E' = Q and then 
applying a purely magnetic formula. The formula of Erber ( 1966| ) with a modification of 
Daugherty & Harding ( 1983| ) correcting its near-threshold inaccuracy is used throughout 
the paper: 

R!{e\B\sine'B) = 

= ci sin S'exp[-C2//(£' sin eisS')], (1) 

where 0'^ = Z(f)', S'), e' is the photon energy; ci and C2 are constant quantities, while / = 
/(e', B') is the near-threshold correction of Daugherty & Harding 1983. Six lines in Fig. ^ 
correspond to six values of the ratio E/B: 0.01, 0.1, 0.2, 0.5, 0.9, and 0.999. Note that the 
free propagation direction deviates from the local B by angle 9''^ which increases with 
increasing contribution of the electric field: 6^^ = arccos(r;^^) = arcsin(£^/i?) ~ E/B. 
Moreover, R increases (A,nfp decreases) monotonically with increasing angle /(?), 77^^). As 
long as E <ti B, magnetic photon absorption rate R remains approximately symmetric 
around the free propagation direction r)^^. 

The directions of t)^^ at various points within the magnetosphere are shown in Fig. |l] 
in the OF (with solid arrows). At the points of photon emission the direction 77^^ overlaps 
with the photon direction 77 since this is where the angle 9b between the photon direction 
and the magnetic field line equals 9b — E/B and consequently R — (Harding et al. 



1978, Zheng et al. 1998). Initially, therefore, the rate R is the same for photons emitted 



from two opposite sites of the polar cap rim which in this picture give rise to the leading 
peak and the trailing peak. But as the photons propagate outward their f} begin to 
deviate from local free propagation directions 77^^ opening thus a possibility for magnetic 
photon absorption and electron-positron pair creation. Three effects are responsible for 
it to happen: (1) curvature of the magnetic field lines, (2) increase of E/B with altitude, 
and (3) slippage of magnetic field lines under the photon's path. For a static-like dipole 
(assumed in section 2) the curvature of lines is symmetric with respect to the magnetic 
axis. Consequently, symmetry is expected between the rates of magnetic absorption for 
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Fig. 2. Directional dependence of the mean free path Amfp for magnetic absorption of a 
photon emitted at the center of the reference frame (with the coordinates [0, 0]) is shown 
for six values oi E/B (denoted with six different line types). The electric field vector 
E is pointing at the right angle to the page, and toward the reader; the E x B vector 
is parallel to the horizontal axis, and pointing to the left. The mean free path Amfp is 
symmetric with respect to the horizontal axis. 

LP photons and TP photons. Any asymmetry between the LP and the TP may occur via 
the effects (2) and (3) only. Whenever the curvature of magnetic field lines happens to be 
relatively large (which is the case at the polar cap rims of pulsars with P ^ 0.1 s) the effect 
(1) significantly dominates over (2) and (3), which means that the resulting asymmetry 
is quite subtle. However, for a relatively small curvature of magnetic field lines (and 
appreciable local corotation velocities /3) the resulting asymmetry becomes pronounced 
[note that these 'favourable' conditions are fulfilled within outer-gap accelerators]. 

We find that the effect (2) ~ the presence of electric field - is crucial for generating 
the asymmetry. The way in which this effect operates can be most clearly assessed by 
inspecting Fig. ^ Let us consider a photon propagation vector r) anchored at the origin of 
frame in the figure. For photons emitted close to the star, both in the LP and the TP, 1] 
initially points along the t)^^ direction which differs only slightly from the direction of B 
in Fig. 1^ since E B. As the LP photon moves away from emission point, its propagation 
vector rotates clockwise in Fig. |^ which reflects the fact that B diverges from ?} due to 
the magnetic field line curvature. At the same time, however, the directional pattern of 
Amfp rotates counterclockwise due to increase in E/B which additionally enhances the 
absorption rate. In the case of the TP photon, however, both its propagation vector fj 
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and the directional pattern of A^fp rotate in the same direction (counterclockwise) with 
respect to B in Fig. || so that the absorption rate is weakened. Thus, for photons within 
the LP the effects (1) and (2) cummulate, whereas for the TP they effectively tend to 
cancel out each other, and the expected asymmetry between the peaks is due to stronger 
absorption suffered by photons within the LP than by photons of the same energy within 
the TP. Accordingly, the high energy cutoff in the LP spectrum will occur at a slightly 
lower energy than the cutoff in the TP spectrum. The slippage (3) affects this picture 
in the way which depends on both the rotation period and the photon position within 
the magnetosphere, but an overall picture remains unchanged. For most rotation periods 
(P > a few ms), the slippage reduces the asymmetry only marginally. For the fastest 
rotators {P ^ 1.5 ms) it enhances the asymmetry by making photons of the TP to 
propagate along the free propagation direction (or, equivalently, along the magnetic field 
lines in the CF). The latter case is presented in Fig. |] where the photon propagation 
direction t) as seen in the OF (dotted lines) and the local free propagation direction 77^^ 
in the OF (solid arrows) are shown for several positions along the photon trajectory in 
the CF. Strong absorption within the LP is anticipated, whereas within the TP fj nearly 
coincides with f)^^. 

Another way to understand this asymmetry is to follow photon trajectories in a 
reference frame (with primed quantities) where the condition 

E' = (2) 

is fulfilled. In such a frame, an asymmetry in R' (cf. eq.(|l|)) for the LP and TP photons 
arises from transformation properties of 6'^ (aberration), e' (red- or blue-shift), and B' . 
One of the reference frames satisfying the condition (|^) is a reference frame of local ExB 
drift. Denoting dimensionless drift velocity (3^ —Ex B/B^, the transformations read: 

Id 

= {^^^) ■ 

In the equatorial plane of orthogonal rotator r^j, = so that rj^ = ^ sin 6b , where the 
signs 'minus' and 'plus' correspond to the leading and the trailing peak, respectively. The 
transformation of propagation angle reduces then to sin6'^ = \rix — f3o\/{l — tIxPd)- The 
Taylor expansion of sin 6*^, e', and B' around /J^ = reveals that a difference between the 
magnetic absorption rates in the locally drifting frame and in the OF results primarily 
from the aberration of photon direction, whereas the changes in e and B are second 
order effects. Obviously, the aberration is asymmetric for the LP and the TP (rjx < 
and r]x > 0, respectively). Fig. |l| presents this 'aberration effect' in the rigidly corotating 
frame, where = is assumed. Photon trajectories in this frame (dashed curves) 
indicate clearly that photons of the leading peak encounter larger B'j^ than photons of 
the trailing peak of the same energy. 
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4. Numerical modelling of the gamma-ray data 



We performed Monte Carlo simulations of curvature-radiation-induced electromagnetic 
cascades developing above a polar cap. The cascade development due to magnetic photon 
absorption accompanied by e^-pair creation and then synchrotron emission was followed 
in a 3D space in order to analyse pulse properties. As an example we choose a model 
with basic parameters of the Vela pulsar: Bpc = 6.8 TG, P = 0.0893 s. In order to meet 
observational restrictions for the Vela, both spectral and temporal, the following general 
requirements within polar-cap scenarios had to be satisfied: 1) a polar-cap accelerator 



should be placed a few stellar radii above pulsar's surface (Dyks et al. |200l[ ); 2) an 



inclination angle a of the magnetic dipole with respect to the spin axis must not be 



large, and the pulsar has to be a nearly-aligned rotator (Daugherty & Harding 1994) 



Recently Harding & Muslimov (1998) proposed a physical mechanism for lifting the polar 
cap accelerator up to 0.5 — 1 i?NS above the surface. However, this altitude is still too 
low to explain the 10 GeV radiation emerging the Vela magnetosphere unattcnuated. 
Therefore, we placed the polar-cap accelerator at the altitude oi Hq — i i?Ns to ensure 
that the magnetosphere is not entirely opaque to curvature photons of energy ^ 10 



GeV (see Dyks et al. 2001 for the detailed model spectral fitting for the Vela pulsar). 



Similarly, Miyazaki & Takahara (1997) achieved the best agreement between the observed 



and their modelled pulse profiles of the Crab pulsar placing the accelerator at /ig = 4 R„ 



To reproduce the observed peak-to-peak separation Ap^^"^ ~ 0.42 (Kanbach et al. 1994) 



we assumed (after Dyks & Rudak 200C) for the inclination angle a and the observer's 



angle Cobs (an angle between the line-of-sight and the spin axis) that a = Cobs = 7.6°. 

Our numerical results are presented in Fig. |^ (a -I- b) . The three columns of Fig. ^ show 
(from left to right): 1) mapping onto the parameter space Cobs vs. cf) of outgoing photons 
with energy e > eumit (where 4> denotes a phase of rotation), 2) double-peak pulse profile 
due to these photons when Cobs — 7.6°, and 3) phase-integrated energy spectrum of these 
photons, with the position of eumit indicated with dotted vertical line. The eight rows 
correspond to 8 different values of £m„it: 1, 10, 10^, 10^ 4•10^ 6•10^ 8•10^ and 10'' MeV 
(top to bottom) . An asymmetry in the double-peak profiles is noticable even though the 
rotator is nearly aligned: at the highest energies, above ~ 6 GeV, the leading peak LP 
is less intense than the trailing peak TP (three lowermost panels in the middle column 
in Fig. |3|b). This is a direct result of stronger magnetic absorption of the LP photons 
comparing to the TP photons. The distribution of these photons in the corresponding 
panels of Cobs vs. <j> (the left column) shows that at viewing angles Cobs larger than 7.6° 
(not allowed due to the fixed peak-to-peak separation of 0.42) the asymmetry in pulse 
profile would be even stronger. This demonstrates an increasing role of rotational effects 
as the distance from the spin axis increases. 
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Fig. 3a. Directional and spectral gamma-ray characteristics calculated for the Vela pulsar 
with the angle a between the spin axis and the magnetic axis set a — 7.6° (nearly aligned 
rotator). Eight rows are shown (the continuation in Fig. 3b), with three panels each. Left 
column shows the outgoing photons of energy e > eumit which are mapped onto the 
parameter space Cobs vs. 0, where ^obs is the viewing angle (between the spin axis and the 
1.0. s) and cf) denotes the phase of rotation. Middle column shows the double-peak pulse 
profile formed with these photons when Cobs — 7.6° is chosen (yielding the peak-to-peak 
separation equal 0.42). Right column shows the phase-averaged energy spectrum (the 
flux level eiv in arbitrary units) for Cobs — 7.6° i.e. the same for all rows. Dotted vertical 
line indicates the part of the spectrum (e > eiimit) which contributes to the corresponding 
pulse profile on the left. The eight rows correspond to 8 consequtive values of eumit: 1, 
10, 10^, 10^ 4 • 10^ 6 • 10^, 8 • 10^, and 10^ MeV (these values are displayed in the panels 
of the right column). 



In the course of magnetic absorption high-energy curvature photons are converted into 
electron-positron pairs which in turn emit low-energy synchrotron photons. Asymmetry 
in the absorption rate as discussed above means, therefore, an identical asymmetry in the 
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Fig. 3b. Continuation of Fig. |3|a. 



pair production rate. Consequently, higher number of low-energy synchrotron photons 
emerges at the LP than at the TP. This is the reason for a dominance of the LP over 
the TP below ^ 100 MeV, noticable in Fig. ||a. Combining the results from both energy 
domains, a characteristic inversion in the relative strentgh of the LP and the TP occurs 
across the gamma-ray energy space. A qualitatively similar inversion of peak intensities 



takes place in the gamma-ray double-pulse of the Vela pulsar (Thompson 2001). 

The beam of synchrotron radiation in our cascades occupies a very narrow range of 
magnetic colatitudes; in other words - it is highly anisotropic. The reasons for this in- 
clude a very limited range of altitudes at which the pairs are created and the effects 
of relativistic beaming. By comparison, curvature radiation below ^ 100 MeV is much 
less anisotropic. Therefore, the prominent peaks visible at e < 100 MeV (two uppermost 
panels of Fig. ^a) consist almost entirely of synchrotron radiation (SR) photons, whereas 
the apparently flat wings outside the peaks (i.e. within the "offpulse" region correspond- 
ing to high altitudes) are composed of curvature radiation (CR) photons. A close-up view 
of the double-peak pulse profile for e > 1 MeV shown in Fig. reveals that the CR wings 
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Fig. 4. Close-up view of the pulse profile for e > 1 MeV from Fig. ga. The curvature- 
radiation wings (|(/)| > 0.3) accompanying both peaks can be now acknowledged. Their 
shape is modelled analytically (shown with thick solid lines). 

are not flat - in fact their intensity decreases with increasing phase |0|; moreover, their 
shapes can be reproduced with analytical means: Spectral power of curvature radiation 
^^g^ well below a characteristic photon energy Scrit oc ^ does not depend on the energy 
of radiating particles 7 but on the curvature radius pcr of magnetic field lines solely. Since 
primary electrons reside within a pulsar magnetosphere for a limited period of time Scrit 
has a lower limit which equals roug hly <, 100 MeV (sec Rudak & Dyks ( |1999D for details). 
Therefore, the wings in the pulse profiles below 100 MeV fall off due exclusively to an 
increase in the curvature radius pcr of magnetic field lines: this proceeds according to the 
following relation 

dPcr 



de 



cx p< 



-2/3 



(4) 



which then leads to the smooth solid lines in Fig. |^. As the photon energy increases and 
the strength of the synchrotron peaks decreases the curvature wings become more and 
more pronounced. They are most noticeable near 100 MeV. Above 100 MeV the wings 
gradually disappear (see the middle column in Fig. ^b) because radiating electrons are 
not energetic enough at high altitudes. At the same time the peaks become narrower - 



an effect noticed by Kanbach et al. (1994) in the EGRET data for the Vela pulsar 



As noted by Daugherty & Harding (1996) the wings within the offpulse region must 
not be too strong within the entire energy range of EGRET if the theoretical pulse 
profiles are to resemble those of the Vela pulsar. We find that the intensity of wings 
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Fig. 5. Pulse profile integrated for e > 100 MeV (left panel) and the phase averaged 
energy spectrum (right panel) are shown for the case of initial energy of primary electrons 
Eq = 2 - 10^ MeV. Note the relatively weaker wings outside the peaks and much softer 
spectrum in comparison with the previous case of Eq = 10^ MeV (Fig. ^ja, third row 
from the top). The data for the Vela pulsar are laid over the model spectrum: boxes 
are the data from the first COMPTEL source catalogue (see Table 3 in Schonfelder et 



al. 2000), and dots (plus an upper limit above 10 GeV) are the EGRET data (Thompson 



et al. 1997). The flux level eF^ is in MeV cm"^ s"^ units. 



relative to the intensity of peaks depends sensitively on the richness of the cascades, i.e. 
on the multiplicity A4 (the number of created pairs per primary electron). The results 
discussed above and presented in Fig. |^ had been obtained for the initial energy of primary 
electrons Eq — 10^ MeV which yielded A4 = 73. By increasing the initial energy Eq up 
to 2 • 10^ MeV the multiplicity reaches M = 830 and the corresponding pulse profile 
at 100 MeV (left panel of Fig. ||) changes notably with respect to its counterpart of 
Fig. I^a. It reveals now a much lower level of wings outside the peaks. Equally important 
is the change in the shape of the phase-averaged energy spectrum which becomes much 
softer by gaining more power in the low-energy range (right panel in Fig. |^). Both new 
features are in rough agreement with the data for the Vela pulsar, contrary to the case 
with Eq — 10*" MeV. It is interesting to note that the association of the broad peaks at 
100 MeV with the relatively hard spectrum (Fig. ^a) on one hand, and of the narrow 
peaks with the soft spectrum (Fig. |^) on the other hand do resemble qualitatively the 
observed characteristics of Geminga and the Vela pulsar, respectively. 

We may now test our model of the double-peak asymmetry by confronting the numer- 
ical results obtained for specific pulsar parameters with the data for real objects. Since 
the effect is induced by magnetic absorption the expected weakening of the leading peak 
with respect to the trailing peak occurs only in the vicinity of the high-energy spectral 
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Fig. 6. The ratio of the number of photons in the traihng peak and the leading peak 
P2/P1 inferred from the EGRET data for Vela is shown in function of photon energy e. 
The results of model calculations performed for three different values of the altitude ho of 
the accelerator are indicated with three lines: dashed, dotted, and solid - for ho — 2i?Ns, 
3i?NS, and 4i?NSi respectively. 



cutoff. Therefore, it is essential to have a good photon statistics also at the highest en- 
ergy bins, i.e. above 1 GeV. As far as the EGRET data are concerned this requirement 
is barely satisfied even for Vela. With these limitation in mind, we consider Vela as the 
only appropriate case to provide the test. We used the EGRET data for Vela to calculate 
the ratio (denoted as P2/P1) of the photon comits in the LP and the TP (denoted as 
PI and P2, respectively). For each energy bin (the energy bins cover the range between 
^ 30 MeV and ~ 10 GeV) we calculated PI (P2) by summing up all photons within 
the range 0ip ± 0.05 (0tp i 0.05) in phase, where (fiip (0tp) is the phase of maximum 
in the LP (the TP) at 100 MeV. Fig. || shows the observational points as well as their 
estimated errors^] along with the results of model calculations performed for three dif- 
ferent altitudes: ho = 2i?NSj 3i?NS, and 4_Rns- The overall qualitative and quantitative 
behaviour of P2 / PI for the EGRET data is very similar to the dependence presented by 
Kanbach et al. 1980| for the COS-B data. The data points certainly can acommodate our 
model. However, to answer the question of whether it would be inevitable to invoke any 
additional processes to reproduce the increase in P2 / PI inferred from the data requires 
better photon statistics at the spectral cut-off and careful statistical analysis. 



^ Detailed analysis of the data will be presented elsewhere (Wozna et al. 2002, in preparation). 
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5. Rotational asymmetry as a function of pulsar parameters 

The asymmetry effects are marginal for nearly aligned pulsars with periods P ~ OT s. For 
the above-described model of Vela, they are noticeable only because of the high altitude 
of the accelerator {ho — 4i?Ns)- However, for highly inclined (a ^ 45°) and fast pulsars 
{P ^ 10 ms) the rotational effects result in asymmetry of considerable magnitude. If 
detected by GLAST, GeV emission from such objects would provide powerful diagnostics 
of the polar cap model. Below we present the magnitude of such asymmetry predicted 
for a wide range of parameters for fast pulsars. 

As a measure of rotational effects we consider the escape energy Eosc which is defined 
as a maximum energy of a photon (and the photon is emitted tangentially to its 'parent' 
magnetic field line; a footpoint of this line has magnetic colatitude 6) for which the 
magnetosphere is still transparent, i.e. for which the optical depth integrated along the 
photon trajectory is less than 1. This energy was calculated with our numerical code and 
the results are shown below. For the sake of comparison with the case of no rotational 
effects we recall a simplified, yet quite accurate^ for magnetic fields weaker than ^ 0.1i?cr, 



analytic formula which (after e.g. Bulik et al. 2000) reads 

£csc{0) ~ • ecsc(^pc), (5) 
where 

lG.v(^)"(jA^)" 

Here is the radial coordinate of the emission point (with spherical coordinates (r„i, 
dm, 4>m) in the right-handed frame with z-axis along the dipole axis) and <C 1 rad is 
the magnetic colatitude of a footpoint of the parent magnetic field line at the neutron 
star surface. For field lines originating at the polar-cap rim one should take then 6 = 6pc, 
where ^^c — arcsin {\f{RNs/R\c)) ~ 1.45- lO^V/P radians, and i?ic = cP/2 TT is the light 
cylinder. Note, that the escape energy of eq.(^ is symmetrical with respect to the dipole 
axis and goes to infinity at the magnetic pole (9 = 0). Such a behaviour is not the case, 
however, when the rotational effects discussed in section 3 are taken into account. 

In Figs. 7, 8, and 9 we present the values of £osc obtained numerically for fast rotators 
with Bpc = 1 TG, and emission points located on the polar cap surface, i.e. rem = -Rns 
was assumed everywhere. Let us begin with the case of orthogonal rotators (i.e. a = 90°) 
- these are shown in Figs. 7 and 8. Here we consider emission points lying along the 
crossection of polar cap surface with the equatorial plane of rotation (hence, the results 
would be relevant for observers located at Cobs — 90°). Fig. 7 shows how £osc varies 
with location of the emission point across the polar cap. The location of each point 



^ For presentation of detailed numerical results see e.g. Harding et al. 



1997 
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Fig. 7. The escape energy £osc of photons from the polar cap surface of an orthogonal 
rotator with Bpc — 10^^ G is shown as a function of normalized magnetic colatitude 
9/9pc of the emission points. The points are assumed to lay along the cross-section of 
the polar cap surface with the equatorial plane of rotation, thus location of each point 
is determined by O/0pc in the range [0, 1], and the magnetic azimuth 0m equal either to 
7r/2 (for the leading half of the polar cap) or — 7r/2 (for the trailing half). Three solid 
lines are labelled with the corresponding spin periods P of 0.1 s, 10 ms, and 1.5 ms. 
Moreover, to illustrate the difference with the case when rotational effects are ignored, 
each solid line is accompanied by a dashed line drawn according to eq.(^. 

is determined by the normalized magnetic colatitude 0/9pc in the range [0, 1] and the 
magnetic azimuth 0,„ equal either to 7r/2 (for the leading half of the polar cap) or — 7r/2 
(for the trailing half). The spin periods of 0.1 s, 10 ms, and 1.5 ms were considered. 

In order to quantify the asymmetry in Ecsc between the leading and the trailing parts 
of the polar cap we introduce the following parameter: 

7^esc(e) = ^, (7) 
£csc 

where 

eJP^ = ecsc (0m = vr/2, 6), 

£*P, = ecsc(0m = -V2, 0). (8) 

For example, 7?.osc(^pc) gives the asymmetry between two opposite points located on the 
outer rim of the polar cap. For spin periods around 0.1 s, typical for the known gamma- 
ray pulsars, 7?.csc(6'pc) remains close to unity; e.g. for the case of P = 0.1 s in Fig. 7, 
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7?.esc(^pc) — 1-3. However, for P smaller than 0.01 s the parameter TZcsd^^pc) becomes 
definitely larger: ^ 2.5 for P = 0.01 s and ^ 30 for P = 1.5 ms. Thus, for pulsars with 
P Si 0.01 s and large inclination angles a we predict a notable difference (in excess of 
half a decade in photon energy) between the positions of high-energy spectral cutoff in 
the gamma-ray emission from the leading and the trailing part of the polar cap. 

Another interesting implication of fast rotation is that ecsc bas finite values for any 
9, including 6 = (the magnetic pole), in contrast to the 'static' case of eq. (|). This can 
understood in the following way: consider emission points with decreasing colatitudc 6 in 
the leading part of the polar cap. As we approach the dipole axis (0 — > 0), Ecsc increases 
because the decreasing curvature of magnetic field lines leads to smaller angles between 
B' and the photon propagation direction r}' in the corotating frame CF. However, photon 
trajectory bends backwards in the CF (see dashed lines in Fig. 1), which implies that also 
photons emitted at 6* = along the straight dipolar axis will quickly encounter 7^ 0, 
thus being subject to the magnetic absorption. Entering now the trailing part of the polar 
cap leads to further increase of Ecsc- This is because magnetic field lines start to bend in 
the same direction as the photon trajectory in the CF (in other words - the efficiency of 
absorption decreases for emission points in the trailing part of the polar cap). Eventually, 
at some point (we denote it as ^o) tbe escape energy reaches a maximum. This is the 
point where the magnetic field slippage along with the aberration of photon direction 
ensure small angles between B' and f/ over large distances in the photon trajectory. 
Therefore, the faster is the rotation, the larger is the colatitude of that point. For 
example, 6'0/6'pc 0.15 for P = 0.1 s whereas 6'0/6'pc ^ 0.42 for P = 0.01 s (see Fig. 7). 
For P = 1.5 ms this maximum occurs ^ close to the outer rim (in the trailing part) and 
therefore huge asymmetry with respect to the leading rim is predicted: 7?.osc(^pc) — 30. 
With further increase of 9 (towards the trailing rim of the polar cap) local magnetic field 
lines start to bend stronger than the photon trajectory in the CF, and this is why £esc 



^ It is easy to reproduce the behaviour of 9o /6pc as a function of P in Fig. 7 with analytic 
formula: The maximal value of B±{r)/ Bpc encountered in a static dipolar field by a photon 



emitted along the local field line at 9o is approximately equal to 0.1 &o (Sturrock 



1971 



also 



Fig.l in Rudak & Ritter 1994 ). This occurs always at ro = 4/3 Rns, regardless the value of 
9o. The angle between the photon propagation direction and the local field line at ro is 
Tpo ~ B±{ro)/B{ro) and therefore ipo ~ 0.1610(4/3)^. To minimize B'j_{ro) as much as possible 
in the case of rotation, the aberration angle due to local linear velocity /3(ro) = 4/3iiNs/^ic 
should be close to ipo- Since /3(ro) <C 1, the aberration angle is ~ l3{ro)- Therefore, we obtain 
the condition 0.1 (4/3)^ 9o ~ 9^^ which gives 

|^«0.08P-i/^ (9) 

Ppc 



This formula overestimates the locations of maxima in Fig. 7 by a factor 2 only. 
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Fig. 8. Escape energy at fixed magnetic colatitude as a function of rotation period. 
The dashed, sohd, and dotted Unes correspond to 6* 2.6°, 8.3°, and 22°, respectively. 
The thin hnes correspond to points on the leading side of the dipole axis (with spherical 
coordinates (i?NS, '^|'^^ (^))- The thick lines correspond to points on the trailing side 
(with (i?NS! ^))- The curves are for emission in the equatorial plane of orthogonal 

rotator. 



should now decrease. However, this decrease never compensates the asymmetry in Scsc 
at 6pc with respect to 9 — 0, i.e. one always ends up with 7?,csc(^pc) > 1- 

Figure 8 presents escape energy e^Fgc and el^^ in function of spin period P. This 
energy was calculated for a fixed position Of^d of emission point, in order to highlight 
its dependence on rotation. We chose three pairs of oppositely located emission points 
at: 9{xd = 2°. 6, 8°. 3, and 22°, which corespond to 6pc for P = 0.1 s, 10 ms, and 1.5 ms, 
respectively. As in Fig. 7, the emission points were placed at the neutron star surface, 
in the equatorial plane of orthogonal rotator. In the case of slow rotation (P ^ 1 s), the 
values of eJPcl^fxd) for the leading point and eosc(^fxd) for the traifing point are practically 
identical, and well approximated by eq.(H). As rotation becomes faster (P around ^ 0.1 
s) £cPc(^fxd) and e|;P(,(6'fxd) start to diverge due to the asymmetric infiuence of E. At 
even shorter periods, below ^ 0.03 s, the maxima in e^^^ are reached (at the values of 
P, which can be reproduced by solving eq.(^ with 0fxd in place of ^o) because slippage 
of field lines starts to be important (see previous paragraph). For increasing 0fxd (2°. 6, 
8°. 3, and 22° in Fig. 8), the asymmetry parameter 7?,osc(6'fxd) decreases for slow rotators 
(P ^ 0.1 s), and it increases for fast (millisecond) rotators. 
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log P [s] 

Fig. 9. Escape energy for the leading and the trailing emission points at the polar cap rim 
as a function of rotation period. Six uppermost lines (thick) correspond to the trailing 
points. Thin lines are for the leading points (the lines for a — 70°, 80°, and 90° overlap). 

In Fig. 9 we present ecsc(^pc) and ecsc(^pc) as a function of spin period P for a set of 
dipole inclinations a. Note, that unlike in Fig. 8, the emission points are now placed at 
9 = 9pc, i.e. at the rim of polar cap corresponding to P. In the case of small inclinations 
(eg. a — 10°), the resulting ratio 7?.csc(^pc) remains close to unity even for millisecond 
periods. The general increase in Ecsc with P increasing, noticeable in Fig. 9, reflects the 
approach of emission points to the dipole axis, where the curvature of magnetic field 
lines is small. The trend is well described by £osc oc \/P as given by eq. (^. In the range 
of spin periods below ~ 0.1 s and with a ^ 45°, the difference between e|P^ and e^^^ 
becomes pronounced, especially for highly inclined millisecond pulsars. If detected by 
GLAST, high-energy emission from such objects would provide an ideal test of the polar 
cap model. 

However, for the asymmetry effects to be detectable, an additional condition (apart 
from short period and large inclination) must be fulfilled: the viewing geometry must 
be of the 'on-beam' type, i.e. the observer's line of sight must cross the narrow beam of 
radiation at the high-energy cutoff - where the magnetic absorption operates. If this is 
the case, the asymmetry in the absorption may be noticeable even in the phase-averaged 
spectra. As an example, we show in Fig. 10 the phase-averaged spectrum calculated for 
a millisecond pulsar with P = 2.3 ms, Bpc = 10^ G, a = 60°, and for Cobs = 60°. For 
this rotator we obtain e[P^ ~ 10'°^ MeV and e^^^ ~ 5 • 10^ MeV at the rim of the polar 
cap. As a result of this rotationally induced asymmetry in magnetic absorption for the 



J. Dyks & B. Rudak: Asymmetry in the lightcurves of EGRET pulsars 19 



T — ' — ' — ' — ' — I — ' — ' — ' — ' — r 




_. . . . I . . . . I . . . . I . . . ._ 

4.0 4.5 5.0 5.5 6.0 
log e [MeV] 



Fig. 10. Theoretical "on-beam" spectrum of high-energy emission from a millisecond 
pulsar with P — 2.3 ms, Bpc = 10^ G, and Cobs — a — 60°. The density distribution 
of primary electrons over the polar cap, which was assumed in this calculation, had a 
maximum at the polar cap rim and was uniform in magnetic azimuth 0m . The magnetic 
colatitude profile of this distribution had a gaussian shape centered at 9 = 6pc, with a 
width a = 0.05 9pc- The initial energy of primary electrons was equal to 2 • 10^ MeV. 
Note the step-like decline near ~ 100 GeV. 



leading and the trailing peak the spectrum at its high-energy cutoff assumes a step-like 
shape: below ~ 10^ MeV the spectrum consists of photons from both the leading and 
the trailing peak, whereas between ^ 10^ MeV and ~ 5 • 10^ MeV only photons of the 
trailing peak contribute to the spectrum; at e ~ e|Pj. the level of the spectrum drops by 
a factor of ^ 2. In these particular calculations of the spectrum, we assumed that the 
density distribution of primary electrons over the polar cap is dominated by an outer-rim 
component (see the captions to Fig. 10 for details of the distribution). 

If, however, an emission from the inner part of the polar cap were to contribute 
considerably to the outer-rim emission, the step-like shape shown in Fig. 10 would be 
smoothed out, because of contribution of many spectra with different values of £osc- In 
such a case the polar cap origin of the observed radiation could be easily revealed by 
noting strong differences between the high-energy spectral cutoffs in different ranges of 
the rotational phase (i.e. phase-resolved spectra would have to be obtained). The first 
obvious candidate to check for this effect (e.g. with GLAST) seems to be J0218-I-4232 - 



the only gamma-ray pulsar among all millisecond pulsars (Kuiper et al. 2000). However, 



this pulsar appears to be a candidate for an 'off-beam' case (see Dyks & Rudak 2002 for 
details). 
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The strength of magnetic field Bp^ practically does not affect the shapes of curves 
shown in Figs. 7-9. The magnetic field only acts as a scaling factor: £osc oc B~^; 
cf. eq. d). 



6. Discussion 

We have shown that pulsar rotation induces an asymmetry in the magnetic absorption 
rate with respect to the magnetic dipole axis. Its consequences are potentially interesting 
in constraining the phase-space of parameters in the polar cap models of high-energy ra- 
diation, provided that very high quality gamma-ray data (e.g. as expected from GLAST) 
are at hand. Its magnitude depends mainly on the linear velocity /3 of the magneto- 
sphere at sites of particle acceleration and magnetic photon absorption. When the region 
of electron acceleration is placed just above the neutron star surface rotation does not 
produce any detectable effects even for relatively fast rotating young gamma-ray pulsars. 
However, it has been argued that at least in the case of the Vela pulsar, such a situation 
is difficult to reconcile with the spectral high-energy cutoff at about 10 GeV (e.g. Dyks 



et al. 2001 ) . We find then that raising the accelerator up to '-^ 4 neutron star radii (in the 



spirit of Harding & Muslimov 1998) above its polar cap produces asymmetric gamma-ray 
pulse profiles even in the case of nearly aligned rotators with a spin period of P ^ 0.1 s. 
The resulting features - softer spectrum of the leading peak and the dominance of the 
trailing peak above ~ 5 GeV - do agree qualitatively with the EGRET data of the bright 
gamma-ray pulsars (Thompson 200 1| ). 



We are far from concluding that the rotation effects alone can account for the observed 
asymmetry in the double peaks of the bright EGRET pulsars. On the contrary - some 
axial asymmetry intrinsic to the region of electron acceleration is inevitable in order to 
explain the double-peak properties at 100 MeV of Geminga and B1706-44, where the 
leading peak is weaker than the trailing peak. Strong deviations of the actual magnetic 



field structure from the pure dipole at the stellar surface (eg. Gil et al. 2002) might be 
responsible for maintaining axial asymmetry at the site of electron acceleration (unlike 
the symmetric initial conditions introduced in Section 2). This in turn would lead to 
electromagnetic cascades whose properties vary with magnetic azimuth. It is important, 
however, that the propagation effects due to rotation work in the right direction, i.e. they 
explain qualitatively the observed weakening of the leading peak with respect to the 
trailing peak. We emphasize that this weakening occurs only in the vicinity of the (phase- 
averaged) high-energy spectral cutoff, where the flux level decreases significantly. 

Another consequence of the magnetic absorption of high energy photons is a noticeable 
change in the separation A^''^'^ between the two peaks in the pulse, taking place near the 



high-energy spectral cutoff (Dyks & Rudak 2000 ). In the model discussed above, with 



electrons ejected only from a rim of the polar cap, the higher energy of photons requires 
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higher emission altitudes to avoid absorption. Therefore, a shght increase in Ap°^'' is 
visible in the three lowermost pulse profiles in Fig. |^b. However, if the emission from the 
interior of the polar cap were included, just the opposite behaviour would occur: A?''^'^ 
would decrease near the high-energy cutoff in the spectrum. This is because in this case 
of a "filled polar cap tube" , the highest energy non-absorbed photons are emitted closer 



to the magnetic dipole axis (see Fig. 2 in Dyks & Rudak 2000). The latter case agrees 
qualitatively with the marginal decrease in peak separation found in the EGRET data 



for Vela (Kanbach |1999D . 

Stimulated by high-quality observations of gamma-ray pulsars anticipated with 
GLAST we analysed in Sect. 5 the importance of rotation-driven asymmetry in mag- 
netic absorption for a broad range of pulsar parameters. A decline in gamma-ray flux 
at high-energy spectral cutoff should inevitably be accompanied by strong changes in 
pulse profiles: whereas at lower photon energies the profile is determined by the density 
distribution of primary electrons over the polar cap and the efficiency of photon emission 
mechanism, in the vicinity of the cutoff it becomes additionally constrained by likely 
high values of the asymmetry parameter 7?.csc(^pc) - the situation anticipated for fastly 
rotating {P < 0.01 s), and highly inclined {a ^ 45°) pulsars. 
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